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The morphology transition from the thermodynamically favorable to the unfavorable phase during 
growth of free-standing copper nanoclusters is studied by molecular dynamics simulations. We give 
a detailed description of the kinetics and thermodynamics of the process. A universal mechanism of 
a solid-solid transition, from icosahedral to decahedral morphology in nanoclusters, is proposed. We 
show that a formation of distorted NC during the growth process with islands of incoming atoms 
localized in certain parts of the grown particle may shift the energy balance between Ih and Dh 
phases in favour of the latter leading to the morphology transition deep within the thermodynamic 
stability field of the former. The role of diffusion in the morphology transition is revealed. In 
particular, it is shown that fast diffusion should suppress the morphology transition and favour 
homogeneous growth of the nanoclusters. 


The study of nanoclusters (NC) growth has recently 
been in focus of many intense research activities [T], be¬ 
cause of their unique properties and also the fundamen¬ 
tally new physical effects that occur in finite-sized sys¬ 
tems [2;. NCs have found broad applications in catalysis 
[3H5] biomedicine or photovoltaic [SHE]- All the applica¬ 
tions require a precise control of the NC growth process 
and understanding of their properties. One of the most 
significant properties of the NC is their morphology. It 
was shown that the favorable structure for copper nan¬ 
ocluster for N<1000 atoms is the Mackay icosahedron 
(Ih) 0 HQ] , followed by the Marks decahedron (Dh) pT] , 
which corresponds to the minimum energy structure for 
1000<N<30000 atoms in the NC. At larger N the clus¬ 
ters optimal structure is FCC [12]. However, experi¬ 
ments demonstrate a significant amount of energetically 
unfavorable morphologies for different synthesis methods 
mm- It is commonly accepted that it happens due 
to kinetics of the growth process. To study the growth 
kinetics is thus of fundamental importance. Baletto et 
al. [15] have demonstrated a possibility of the morphol¬ 
ogy transition in molecular dynamics simulations of sil¬ 
ver nanoclusters and explained its mechanism PH- Th ey 
claimed that the transition from icosahedral (Ih) to dec¬ 
ahedral (Dh) morphology goes through a melting of the 
cluster and a formation of the new morphology from the 
melted cluster. Lan et al. m proposed a qualitatively 
different mechanism for the transformation from Dh to 
Ih, where morphology goes through a solid-solid transi¬ 
tion without the formation of an amorphous phase. On 
the other hand, it is known that the smaller the cluster 
size is the more favorable the Ih phase should be. This 
means that in a real growth process the transition should 
go from Ih to Dh and it is very unlikely to go backwards. 
Thus it is much more relevant to study the Ih-Dh tran¬ 
sition. 

In this Letter, we report the mechanism of the solid- 
solid morphology transition that we have discovered in 


copper nanoclusters. Also, we provide a detailed ex¬ 
planation of the conditions that induce the morphology 
transition, or on the contrary preserve the layer-by-layer 
growth. Also, the decisive role of diffusion in the mor¬ 
phology transition is exhaustively explained. 

Note that the main interest in MD simulations of NCs 
addressed low-energy nanocluster growth techniques. 
However, recent developments of new NC synthesis meth¬ 
ods, e.g. employing pulsed highly ionized plasma mm, 
have drawn attention to higher energies during the 
growth. The use of a plasma environment has several 
important advantages such as an increased growth rate 
and a wide choice of possible materials to grow. It has 
already been shown that the energy of the growth pro¬ 
cess significantly affects the kinetics of the growth process 
m , a fact that should be accounted for in simulations. 
For instance, it was shown that the Coulomb interaction 
between the NC and impinging ions during the growth 
process will influence the angular distribution of the ve¬ 
locities of the impinging atoms. Consequently the angu¬ 
lar distribution will affect the diffusion and local heating 
of the cluster surface. 

We used the embedded-atom method (EAM) potential 
with Foiles parametrization m to simulate the growth 
of copper NCs. We start our growth simulations from a 
seed. Based on the Baletto’s analysis of structures and 
magic numbers for clusters [T2[, an icosahedral seed con¬ 
sisting of 147 atoms was chosen as our growth seed. In 
the simulation new atoms were randomly generated ev¬ 
ery 100 ps around the cluster on a sphere with radius 
13 A. That is, slightly outside the cutoff radius of the 
potential from the NC. Though the growth rates in our 
simulations are much higher than the experimental ones, 
we employ a thermostat to cool the cluster down to the 
correct temperature before a new particle impinge, this 
provides the correct conditions for the growth. The final 
structure was analyzed after that 300 atoms were added 
to the seed. The temperature of the nanocluster was con- 
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trolled by a Nose-Hoover thermostat m and the speed of 
the impinging atoms corresponded to either the temper¬ 
ature of the cluster or a kinetic energy of 1 eV. These two 
cases were chosen to represent the inert-gas-aggregation 
(IGA) growth process or the conditions of the growth in 
a plasma, respectively m- 

We performed simulations in LAMMPS [23] for dif¬ 
ferent temperatures, incident velocities and angular dis¬ 
tribution of the impinging atoms [20 ]. Every considered 
temperature, velocity and angular distribution of the im¬ 
pinging atoms were studied in series of 50 simulations. 
Analysis of the results showed significant fractions (up to 
30% ) of non-icosahedral clusters among the final clus¬ 
ters. All results are assembled in Table I. Normal inci¬ 
dence and full distribution correspond to particles grown 
from atoms with an energy of 1 eV, simulating the NC 
growth in plasma. We have considered two models of the 
growth, with atoms impinging on the surface with nor¬ 
mal incidence and with the angular distribution derived 
in [20] . Results denoted as IGA correspond to a respec¬ 
tive process simulated with the assumption that the im¬ 
pinging particles have a thermal energy (0.03 eV). Note 
that, as it was mentioned early, the thermodynamically 
favorable morphology for clusters of this size is icosahe- 
dral. Thus the kinetics of the growth process induces the 
morphology transition. However, the majority of the fi¬ 
nal clusters are still icosahedral, which rises the question: 
What is the difference between these two outcomes? 

To address this question, we analyze in details one typ¬ 
ical simulation, where the structure changes from icosa¬ 
hedral to decahedral during the growth. Fig.l shows co¬ 
ordination of atoms calculated with Bond-angle analysis 
(BAA) [23. The method was developed by Ackland and 
Jones [25i to distinguish fee, hep and bee coordination 
structures. From the bond vectors of the central atom 
the histogram of the bond angle cosines is computed and 
then used to determine the structure type by the heuris¬ 
tic decision rules. Initially the cluster consists of 147 
atoms arranged in a perfect icosahedron with all atoms 
HCP-coordinated. However, all surface atoms are identi¬ 
fied as disordered by the BAA method. Thus we see that 
« 50 interior atoms are HCP-coordinated and ^100 sur¬ 
face atoms are disordered. At the size of « 250 atoms 
the cluster undergoes the morphology transition from 
icosahedral to decahedral. At that point we see how the 
number of HCP-coordinated atoms drops down and the 
number of FCC-coordinated atoms abruptly increases. 
It clearly shows that part of the HCP-coordinated atoms 
change their local structure to FCC-coordination, since 
we don’t see any significant change in disordered atoms. 
After the transition FCC-coordinated atoms dominate 
over the HCP-coordinated ones for the remaining of the 
growth simulation. This means that there were no more 
morphology transition. 

In order to understand what induces the transition we 
analyzed the potential energy of the cluster as a function 
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Figure 1: Bond-angle analysis of the cluster structure. The 
figure shows how the number of HCP, FCC and disordered 
atoms changes with time, or equivalently the cluster size, dur¬ 
ing the growth process. 


of time. Fig.2 shows the typical case described in the 
proceeding paragraph where the morphology transition 
takes place between incident atom 101 and 102, i.e. when 
the cluster has 248 atoms. We see that from 20 till 35 
ps, after atom 101 were introduced, the potential energy 
of the cluster fluctuates. This means that the structure 
is far from perfect and that atoms on the surface change 
their positions and thereby change the potential energy 
of the system. However, after the transition, the struc¬ 
ture falls into a potential well, where it minimizes the 
potential energy and thus we see no more fluctuations of 
the energy. It is worth mentioning that the full transi¬ 
tion takes approximately 2 ps and the whole time range 
of Fig.2 corresponds to one point in Fig.l. The energy 
between these two states is 1.8 eV. Also one can see in 
Fig.2 that the surface area of the cluster changes signifi¬ 
cantly at the transition. The fluctuations of the potential 
energy correspond to changes of the surface area until the 
clusters undergoes the morphology transition. Then the 
surface area promptly increases while fluctuations vanish. 

We applied the Nudged Elastic Band (NEB) method 
to obtain the height of the potential energy barrier for 
the transition from icosahedral to decahedral morphol¬ 
ogy, Fig.2a. This allows us to plot a Potential Energy 
Surface (PES) of this particular cooper nanocluster con¬ 
sisting of 248 atoms, Fig.3. The icosahedral phase corre¬ 
sponds to a minimum with a potential energy of -790.5 eV 
and the decahedral phase corresponds to a deeper local 
minimum with a potential energy of -792.3 eV. The bar¬ 
rier the cluster needs to overcome in order to minimize 
its energy is 0.8 eV. Besides, we calculated the barrier 
heights for different sizes of NCs. The cluster consist¬ 
ing of 185 atoms has a 0.6 eV barrier and the 376-atom 
cluster has to overcome a barrier of 2 eV. 
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Figure 2: Data from the same simulation as in Fig.l but only 
for a short time period when the cluster consist of 248 atoms, 
(a) NEB analysis of the transition from the Ih to the Dh 
phase; (b) potential energy of the cluster as a function of time; 
(c) dependence of the cluster surface area on the simulation 
time. 



ICOSAHERDAL DECAHEDRAL 


Figure 3: The potential energy surface for the same 248 atom 
copper nanocluster as in Fig.l and Fig.2 

It is known, that icosahedron is the favorable morphol¬ 
ogy in many noble metals m at small sizes due to its 
small surface area (see Fig.4a). Thus icosahedron mor¬ 
phology minimizes the surface area, but causes relatively 
large internal stress. On the contrary, the decahedron 
morphology has a smaller internal stress but a larger sur¬ 
face area. If the impinging atoms during the growth are 


Table I: The fraction of the clusters that suffered a morphol¬ 
ogy transition for different temperatures. Normal incidence 
(NI) correspond to particles grown from atoms impinging on 
the cluster with normal incidence and a kinetic energy of 1 eV. 
Full distribution (FD) corresponds to particles grown from 
atoms with kinetic energy of 1 eV and a distributions of inci¬ 
dences from normal to grazing. Particles grown from atoms 
with a kinetic energy of 0.03 eV represent the IGA process. 


Temp, K 

NI, % 

FD, % 

IGA,% 

400 

28 

16 

26 

450 

20 

16 

22 

500 

10 

10 

18 

600 

8 

8 

12 

700 

4 

2 
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not able to diffuse far enough to form an uniform new 
layer, and thus instead form islands on the facets of the 
icosahderon Fig.4b, then the surface area significantly in¬ 
creases. Thus, the surface area of the non-perfect icosa¬ 
hedron is large as well as its internal stress. A transi¬ 
tion to the decahedral phase minimizes the internal stress 
though the surface area increases. Still, as one sees in Fig. 
3a, the potential energy of the decahedral phase is lower 
than that of the imperfect icosahedron. 

Another way to minimize the energy is to redistribute the 
impinging atoms more uniformly, thus preventing islands 
to form on the cluster surface and strong expansion of the 
surface area. Statistics of the simulated clusters (Table 
I) witnesses that the higher the temperature the lower 
the probability of the transition to the decahedral phase. 
This fact can be explained by the rate of diffusion, which 
is higher at high temperatures. Moreover we have ob¬ 
served that small islands may diffuse, thus redistributing 
extra atoms more uniformly on the surface of the cluster 
and thereby decreasing the surface area. At the same 
time, it is clear that growth processes that restricts dif¬ 
fusion like growth with normal incidence particles or the 
IGA process, where the diffusion is limited by the low 
energy of the impinging particle, have noticeably higher 
statistics for the morphology transitions. All these ob¬ 
servations demonstrate how important the diffusion is in 
the morphology transition. But exactly how does a non¬ 
perfect icosahedron transform into a decahedron? 

Detailed observations of hundreds of our NC growth 
simulations allowed us to identify a pattern and describe 
the mechanism of the morphology transitions from Ih to 
Dh phase. This mechanism has very simple foundations 
and thus seems very elegant. The basic idea is the fol¬ 
lowing: the formation of islands on the surface disturbs 
the perfect facets of the icosahderon (Fig. 4a). This 
disturbance can be localized in a certain area of the NC, 
which could be considered as a "belt" formed by the extra 
atoms, Fig4b. Above and under this "belt" one can dis¬ 
tinguish two pentagonal pyramids with much less disor- 
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top pyramid 




Figure 4: Mechanism of the morphology transitions from Ih to Dh phase: a) A perfect icosahedron; b) impinging atoms form 
islands on the surface of the cluster, thus increasing the surface energy, the most disturbed facets could be described as located 
in some kind of "belt"; c) the bottom and top pentagonal pyramids twist with respect to each other and a fractional part of the 
islands on the facets incorporates into the cluster layers; d) aligned edges of the new 100 facets and with all atoms in between 
the slices rearrange into a FCC structure. 


der. These pyramids do not require much rearrangement, 
because they should remain the same in the decahedral 
cluster. They only need to be twisted with respect to 
each other in order to align their edges and form the new 
rectangular facets of a decahedron Fig4c-d. During these 
relative rotations the atoms in the "belt" undergo major 
rearrangements. This includes the formation of FCC- 
coordinated atoms from HCP ones and incorporation of 
the extra atoms from the surface into new layers. 

Baletto et al. m suggested that this transition in sil¬ 
ver nanoclusters goes through a melted phase of the clus¬ 
ter and that the new phase is formed from an amorphous 
structure. However, Lan et al. 03 gave very serious ar¬ 
guments in favor of a solid-solid transition. On the other 
hand, Lan et al. showed a transition in "perfect" clusters 
(cluster with magic number of atoms) from the least fa¬ 
vorable to the most favorable phase in a transition that 
has been driven by pure thermodynamics. Thus they did 
not account for the influence of the kinetics, which is be¬ 
lieved to be the morphology-determining process at real¬ 
istic conditions. Koga et al. suggested a mechanism 
similar to ours from experiments of gold nanoclusters 
growth. We have analyzed more than 700 simulations 
and we have shown that the mechanism of the morphol¬ 
ogy transition remains the same for all considered sizes, 
temperatures and even energies of the impinging parti¬ 
cles. This allows us to claim that this is the universal 
mechanism of the solid-solid Ih-Dh morphology transi¬ 
tion for all copper nanoclusters. That gives us a very 
rigid reason to expect it to be the universal mechanism 
for all metal nanoclusters. 

In summary, we have reported the mechanism of solid- 
solid transition form Ih to Dh phase in metal nanoclus¬ 
ters, and suggested an explanation of the cause of a mor¬ 
phology transition in Ih nanoclusters deep within their 
thermodynamic stability field. The origin of the transi¬ 
tion has been associated with a formation of distorted 
NC during the growth process with islands of incoming 


atoms localized in certain parts of the grown cluster. The 
distortions change the energy balance between Ih and Dh 
phases, leading to the changes of the morphology of the 
NC. We revealed the role of diffusion on the morphology 
transition and showed how different facets of the cluster 
will transform during the transition. The fundamental 
understanding of the morphology transition mechanism 
allows us to suggest that varying the diffusion length in 
the growth process can be used to influence the clus¬ 
ter morphology. Considering the fact that the transition 
barrier increases with the increase of the cluster size we 
expect that once the morphology is set at a certain large 
size it is very unlikely to change. A fact that again seems 
very convenient for tailoring the morphology of the grown 
nanoclusters. 
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